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Abstract. We prove the existence and uniqueness of stationary spherically 
symmetric positive solutions for the Schrodinger-Newton model in any space 
dimension d. Our result is based on an analysis of the corresponding system 
of second order differential equations. It turns out that d = 6 is critical for the 
existence of finite energy solutions and the equations for positive spherically 
symmetric solutions reduce to a Lane-Emden equation for all d > 6. Our result 
implies in particular the existence of stationary solutions for two-dimensional 
self-gravitating particles and closes the gap between the variational proofs in 
d = 1 and d = 3. 



1. Introduction 
We consider the Schrodinger-Newton equations 

(1.1) ii/; t + Aip — jVip = 0, AV=\ip\ 2 

on M. d which is equivalent to the nonlinear Schrodinger equation 

(1.2) iip t + M + 7 (G d (W) * M> = 

where Gd ( | cc | ) denotes the Green's function of the Laplacian on R d . Of physical 
interest are solutions having finite energy E and particle number (or charge) iV 
given by 

(1.3) Ety) = \( \^{x,t)\ 2 dx-l [[ G d (\x-y\m Xl t)\ 2 \^(y,t)\ 2 dxdy 
and 

(1.4) 2Vty)= / \^{x,t)\ 2 dx, 
respectively. 

In the present work we study in the attractive case 7 > the existence and 
uniqueness of spherically symmetric quasi stationary solutions of the form 

(1.5) i/j{t,x)=uU\x\)e- lut , uM)>0, Iim u u {\x\) = 0, 

|a:| — >oo 
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which we call ground states. For solutions of the form (|1.5|) we have V(t, x) = 
Vu}{\x\) and u w (r),u w (r) satisfy the following system of ordinary differential equa- 
tions: 

We suppose that u u (0), v u (0) are finite and u^(0) = v' w (0) =0. The latter equation 
implies that v' u > and therefore for solutions u u vanishing at infinity we have 
w — 1 v lu{0) > 0. By rescaling u(r) — Au u (r/a), V(r) — B{v u {r/a) — w/7) + 1 with 

G G 



we obtain the universal equations 

„ d - 1 



(1.7) 



T 



v! = (V- l)u 



subject to the initial conditions 

(1.8) u(0) =u Q e R+ u'(0) = 0, 7(0) = 0, V'(Q) = 0. 

Here d > may be regarded as a continuous parameter. By analyzing the 
solutions of the above initial value problem we shall prove the following result 
about the existence and uniqueness of ground states: 



Theorem 1.1. For any d > the system (|1.7j) subject to the initial conditions 
(|1.8p admits a unique solution (u,V) such that u(r) > 0, u (r) < on ]0, 00 [ and 

(1.9) lim u(r) = 0. 
If d = 6 ? £/ien w(r) zs explicitly given by 

(1.10) u ( r ) = fi + ?_ 

In addition, for all d > 6 t/ie relation u(r) = 1 — V(r) ZioWs and </ie function u(r) 
solves the Lane-Emden equation 

/ 1 1 1 *\ it , d — 1 , 2 

(1.11) it H u = —u 

r 

subject to the initial conditions u(0) = 1, u'(0) = 0. 

The decay properties of these solutions will imply that they have finite energy 
and particle number if and only if d < 6. 

In the physical and mathematical literature the Schrodinger-Newton system in 
three space dimensions has a long standing history. With 7 designating appropriate 
positive coupling constants it appeared first in 1954, then in 1976 and lastly in 1996 
for describing the quantum mechanics of a Polaron at rest by S. J. Pekar pQ, of an 
electron trapped in its own hole by the first author [2] and of self-gravitating matter 
by R. Penrose [3J. In 1977, E.Lieb [5] showed the existence of a unique ground state 
of the form (jl.5p in three space dimensions by solving an appropriate minimization 
problem. This ground state solution u^(x), uj > is a positive spherically symmetric 
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strictly decreasing function. In [4 , P.L. Lions proved the existence of infinitely 
many distinct spherically symmetric solutions and claimed a proof for the existence 
of anisotropic bound states in [5]. 

While Lieb's existence proof can be easily extended to dimensions d — 4 and 
d = 5, the situation is unclear for lower dimensions due to the lack of positivity 
of the Coulomb interaction energy term. For the one-dimensional problem this 
difficulty has been overcome recently in [6] and the existence of a unique ground 
state of the form (jl.5p has been shown by solving a minimization problem. The 
two-dimensional problem, however, remained open and so far only numerical studies 
are available indicating the existence of bound states, see e.g. [7]. Our main result 
proves the existence of such solutions. 

From the variational point of view the critical dimension d = 6 is related to a 
Hardy-Littlewood-Sobolev inequality of the form 



for a positive constant C. Our solution (|1.10|) is indeed an optimizer in this in- 
equality. Instead of proving this inequality directly we deduce it together with 
the optimal constant by simply combining two inequalities of [8] (see Appendix). 
From the ODE point of view the system of ordinary differential equations has a 
conformal invariance for d — 6 which leads to a one-dimensional autonomous sys- 
tem with a Yukawa-type interaction. In addition, the problem of finding a positive 
solution can be reduced to solving a (conformally invariant) Lane-Emden equa- 
tion [9] , [10] , [UJ . We do not make use of this property but we believe that this 
observation may be useful for further studies of these equations and we give the 
corresponding autonomous system in the Appendix. 

To prove the main result we use a shooting method which various authors have 
successfully applied to existence and uniqueness of solutions in boundary value 
problems for second order nonlinear differential equations [12], [13], [14], [15] . 
Our paper is organized as follows: In Section 2 we employ a shooting method to 
prove the existence of ground states (theorem 12. 3ft . In Section 3 study their decay 
properties to prove uniqueness by analyzing the Wronskian of solutions (theorem 
I3.6[) . Finally, in Section 4 we prove the final part of our main theorem including 
the explicit solution for d = 6 (theorem 14. \\ . 



We begin our study with the discussion of some general properties of solutions 
of (|1.7|1 with initial values (jl.8p . Standard results will guarantee local existence 
and uniqueness of solutions, their continuous dependance on the initial values as 
well as on the parameter d and their regularity. As a consequence of local existence 
and uniqueness solutions cannot have double zeros. We shall frequently apply these 
properties in the sequel as well as the following integral equations for u' and V: 



(1.12) 




2. Existence of ground states 



(2.1) 
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Viewed as a mechanical system we can associate an energy to (|1.7|) given by 

(2.2) 2£ = u' 2 + u 2 + ^V' 2 - Vu 2 
which satisfies 

£' = u' 2 - — V' 2 . 

r 2r 

Therefore £ is a constant of motion if d = 1. However, we shall not use this property 
in the present work. 

For the initial condition uq > of the solution (u, V) we consider the following 
mutually disjoint sets: 

Definition 2.1. 

(2.3) TV = {u e M+ : 3r > such that u{r ) < and u'(r) < on ] 0, r ] }, 

(2.4) g = {u e M+ : u > 0, lim u(r) = 0}, 

r too 

(2.5) V = {u S M+ : 3n > such that u'(n) > and u(r) > on ] 0, n] }. 

In order to see that Q and P are disjoint note that since u"(0) = —duo all 
solutions start strictly decreasing. Therefore any solution with initial condition uo 
in V has a local minimum before n where V > 1. From (|2.ip we deduce that V 
is strictly increasing. Therefore u' will remain positive and bounded away from 
zero after t\ and u becomes unbounded. If uq ^ AfUP, then u > and u' < 0. 
Therefore it has a limit as r tends to infinity which must be zero. This implies 
AfUGUV = M.+ . From the continuous dependance on initial values we deduce that 
TV, V are open sets. 

Our main result theorcm ll.il states that Q consists of exactly one element. Ob- 
viously, Q is nonempty if N and V are nonempty which we show in the following 
lemma. 

Lemma 2.2. For any d > the sets Af and V are nonempty. In particular, 
]0,1[C/V. 

Proof. We consider the function 

(2.6) <P = u + V-l. 
It satisfies the differential equation 



4>' = u<f) 



and admits the Taylor expansion 



0(r) = K-l)(l + ^)+o(r 2 ). 

Let Uo < 1 and suppose Uq ^ Af. Then <j) is negative and strictly decreasing for all 
r > sufficiently small. By hypothesis u remains strictly positive. Consequently, 
<j> cannot have a critical point since then <j>" = u<f> < which is impossible. We 
conclude that cj>(r) < <fi(0) = uq — 1 for all r > or equivalently, u(r) + V(r) < uq. 
Hence uq e Q and V(r) < uo for all r > 0. Since V is always strictly increasing 
Voo := lim V(r) exists and Voo < uq < 1. We consider the function z :— — 
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Since uo G Q it follows that z exists for all r > and z(r) > for all r > 0. It 
satisfies the differential equation 



z =z 2 - - z+l-V. 

r 

Choose f such that < i/2(l — Voo) for all r > f. Then for r > f we have 

z > \z 2 + (iz 2 - ^/2(l-V 00 )z + l~vjj> X -z 2 . 

This implies that z blows up in finite time which is impossible. Hence uq € TV. 

Next we want to show that uq € V for uq sufficiently large. Suppose on the 
contrary that V is empty and let uo > 1. Denote ]0, Rq[ the maximal interval 
where u > and v! < 0. Therefore from equation (|2.I|) for V' we obtain the 
bounds 

a a 

Integrating these inequality and using again that u is decreasing yields the following 
estimates for V: 

u(r) 2 r 2 , , vlr 2 

By a similar reasoning as before we see that the function <fi defined in l|2.6p is strictly 
increasing on ]0,i?o[- Hence 

u(r) > u - V(r) on]0, R [. 



d-1 



Let ro = y/2d/uo. Inserting the upper bound for V we get 

r 2 

u(r) > u (l 011 ]°) r o[ 

r o 

and ro < Rq- We want to show that u'(ro) > provided uo is sufficiently large 
which yields the desired contradiction. Using our bounds on u and V in (|2.1|) we 
obtain 



1 



''0 



u'M = -ttt / (V(r) - l)u{r)r d - x dr 





ro 



-2d4' l L () d 



> ^{du j\l-* l ) 3 s d+1 ds-1 

We conclude that u'(ro) > for uq sufficiently large which contradicts the assump- 
tion that V is empty. □ 

Hence we have proved by the preceding lemma the existence of ground states: 

Theorem 2.3. For any d > 1, the set Q is nonempty, that is there is a solution 
(u, V) (|1.7|) subject to the initial conditions (|1.8j) such that u(r) > 0, u'(r) < on 
]0, 00 [ and lim u(r) = 0. 
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Proof. It remains to prove u(r) > 0, u'(r) < 0. The first property follows from the 
fact there are no double zeros. If u has a first critical point r% > 0, then V(ri) > 1 
and since V is strictly increasing (see (I2.1P ) it follows again from (|2.1[) that u(r) > 
for all r > r\ which is impossible. Hence v! < 0. □ 

3. Uniqueness of ground states 

In this section we prove that Q has exactly one element. First of all, we show that 
if Q had more than one element the corresponding solutions cannot cross. This is an 
immediate consequence of the following lemma which states that any two solutions 
of the initial value problem (ll.7|) . ()1.8|) cannot cross as long as they stay positive. 

Lemma 3.1. Let 1*2(0) > ui(0) > and suppose that u 2 (r),ui(r) exist on [0,R] 
such that u\(r) > on [0, R). Then U2{r) > u\{r) for all r £ [0, R\. 

Proof. We consider the Wronskian of U\ , u 2 defined by 

(3.1) w(r) = u' 2 {r)u\(r) — u 1 {r)u2{r). 
Then w satisfies the differential equation 

(3.2) w' + — - w = (V 2 - Vx)uiu 2 . 

r 

Suppose there is f G [0, R] such that u 2 (r) > u\{r) on [0, f[ and u\{f) — u 2 (f) > 
0. Then 

w(f) = (u' 2 (f) — u[(f))ui(f) < 0. 
On the other hand we have 

Vi(r) ~ V[{r) = f(ul(s) - u?^))^" 1 ds > 

r Jo 

on ]0, f] and therefore V 2 (r) > V\ (r) on ]0, f]. We conclude then from the differential 
equation (|3.2[) for w that u>r d_1 is strictly increasing on ]0,f] and since w(0) = 
we must have w(f) > which is the desired contradiction. □ 

Remark 3.2. From the no-crossing property stated in lemma 13.11 it follows im- 
mediately that TV, V are intervals. More precisely, M =]0, a[,V =]b,oo[ with 
0<a<6<oo[. Uniqueness of ground states is then equivalent to a = b. 

The important conclusion from lemma 13.11 is that two different ground state 
solutions cannot intersect. From the differential equation (|3.2[) for their Wronskian 
w(r) we see that w{r)r d ~ 1 is a nonnegative strictly increasing function. However, 
we shall prove in the sequel that w{r)r d ~ 1 vanishes at infinity which yields the 
desired contradiction. Therefore we have to analyze the decay properties of ground 
states at infinity. 

Since V is always strictly increasing Voo := lim V(r) exists (including the case 

r > 00 

Uoo = +00) and V{r) < for all r > 0. 

Lemma 3.3. Let uq S Q. Then 1 < uq < Voo- In particular, uq — 1 if and only if 
Voo = 1 and in this case u = 1 — V . 

Proof. Since Af c]0, 1[ by lemma [2~2l the inequality 1 < uq is obvious. To prove the 
second inequality it is sufficient to consider the case < +00. We consider the 
function £ defined by 

£ = U + V -Voo. 
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Obviously, £(0) = uq — Voo, £'(0) = and lim £(r) = 0. The function £ satisfies 

r — >oc 

the differential equation 

£"+ — £' = «£ + (Voo- 
r 

and admits the Taylor expansion 

t( \ u j u o(u - f )r 2 2 
f (r) = u Q - Foo H — V o(r ). 

Suppose uq > Voo. Then £(r) is positive and strictly increasing for r > sufficiently 
small. Therefore £ must have a critical point ri with £(ri) > Uo — V» and < 

which is impossible since u > 0. 

For the initial condition u(0) = 1 we have thanks to the uniqueness of solutions 
for the initial value problem \1.7\ . (| 1 .8(1 that u = 1 — V. □ 



In the following lemma we determine the asymptotic behavior of ground states. 
Lemma 3.4. Let uq G G ■ Then 



lim — = - x/Voo - 1. 



+ oo u 



Moreover, if 1 < < oo, then for any K G] 0, yVoo — 1 [ .■ 

lim sup u(r)e Kr < oo. 

r >oo 

Proof. First off all, let < oo. We consider the function z := — — which is well 
defined for all r > and satisfies the differential equation 

z' = z 2 z + l-V. 

r 

Now choose f such that — - < i-v/V^ for all r > f. Consider the direction field in 



2 

the (r, z) plane for the preceding differential equation. In the set r > f, z > 2y / V oc 
we have 

z' > ^z 2 + (iz 2 - ]^^/V~z + l-Voo)> \z 2 + 1. 

It follows that, should z(r) ever enter this region, it would blow up at finite time 
after r which is impossible. Hence z remains bounded. This also implies 



lim u (f) = 0. 

r 'too 

Therefore we may apply l'Hospital's rule. We obtain 

u" /d—1 \ 
lim z 2 = lim — = lim I z + V — 1 ) = Voo — 1< 

r *oo r >oo u r ^oo \ r / 

Finally, if Voo is infinite, then z is also unbounded since otherwise applying l'Hospital's 
rule as above yields the desired contradiction. This proves the first part of the 
lemma. 

Now let Voo > 1. Then for any k g] 0, VVoo — 1 [ and r sufficiently large, —^>k 
and the proof is completed by integrating this inequality and taking exponentials 
on both sides. □ 
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Remark 3.5. If Voo = +00 the asymptotic behavior of ground states can be given 
more precisely. Indeed, by analyzing the differential equation for Z := — ii-V -1 / 2 
and taking into account that lim V' /V — it can be easily shown by mimicking 

r >oc 

the proof of the preceding lemma that lim Z = 1 . 

r >oo 

Now we are in position to prove our uniqueness result: 
Theorem 3.6. The set Q has exactly one element. 

Proof. Let iti(O), 112(0) £ Q such that 112(0) > iti(O). By lemma [3~T1 the correspond- 
ing solutions u\,u 2 cannot intersect and we have u 2 [r) > u\{r) > for all r > 0. 
From the differential equation (|3.2[) for their Wronskian w(r) we see that w(r)r rf_1 
is a nonnegative strictly increasing function since 

(id/- 1 )' = (V a - Vi)uiu 2 r d - 1 > 

and w(0) = 0. On the other hand, we claim that 

lim wr d ~ x = 0. 

r >oo 

2 2 

Indeed, by lemma [3~31 we have lim V 2 {r) > 1. Trivially, V(r) < U2{ "> r . From the 
integral equation (|2.ip for u\ 

u'^ry- 1 = f r (V 2 (s) - l)u 2 (s) s d - 1 ds 
Jo 

and the decay properties of u 2 given in lemma 13^1 it follows then that it 2 r rf ~ 1 and 
u 2 r d ~ x are uniformly bounded. Therefore 

\w(r)r d - x \ < luilKr^l + KHuar^- 1 ! <ci|ui| + c 2 K| 

for some positive constants ci , c 2 which concludes the proof. □ 



4. Further properties of ground states 

For the initial condition u(0) = 1 we have thanks to the uniqueness of solutions 
for the initial value problem (|1.7p . (|1.8[) that u = 1 — V and therefore we have to 
solve the following initial value problem 

(4.1) u" + — -u' = -u 2 , u(0) = 1, u'(0) = 0. 

r 

This is the d-dimensional Lane-Emden equation. The behavior of its solutions has 
been widely studied in the mathematical literature. However, in the following we 
will give alternative proofs of the results relevant for our work. Analyzing the 
behavior of the solutions of the initial value problem (|4.ip we prove the final part 
of our main result: 

Theorem 4.1. If d > 6, then 1 6 Q and if d < 6, then 1 € J\f. In particular for 
d = 6 the ground state solution is explicitly given by 

(4-2) u(r)= (l + Q 2 ■ 
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Proof. Since ]0, l[c TV by lemma l2~2l and V is open, 1 ^ V for any d > 0. Now let 
d > 6. We introduce a new Lyapunov function L(r) defined by 

(4.3) L(r) = E(r)r d + ~u{r)u' (r)/" 1 . 

o 

Then L(0) = and 

L'M = ^(rV^ 1 <0 

6 

since d > 6. Hence L(r) < for all r > 0. We suppose that 1 € TV. Then there 
exists r > such that u(r ) — 0. Computing L at this point we get 

L(ro) = \u l2 (ro)4 > 

which is impossible. The explicit solution (| 1 . 10[) for d = 6 is readily verified (see 
also Appendix A). 

Let d < 6 and assume 1 E G- We may then use the Milne variables (see 
e.g. m, P) 

— ru' —ru 



u w 

which are well defined for all r > 0. Indeed, y(Q) = 0, z(0) = d and y, z > for all 
r > 0. They satisfy the differential equations 

y' = ^{2-d + y + z), z'= ~(d-2y-z). 
r r 

If d < 2, then the first differential equation implies that y blows up in finite time 

which is impossible. If d > 2, then y and z remain bounded, which implies that for 

all r sufficiently large 

u(r) < Or' 2 , -v! < Cr~ 3 

for an appropriate constant C > 0. We consider now the Lyapunov functional L(r) 
defined in (|4.3[) . Since d < 6 The bounds on u, v! imply that 

lim L{r) = 

r — >oo 

On the other hand, L(0) = and 

L '( r ) = -^—l u a {r)r i - 1 >0 
6 

which yields the desired contradiction. □ 
Corollary 4.2. Let d < 6 and uq G G ■ Then Voo > uq > 1. 

Proof. The inequality uq > 1 is an immediate consequence of the previous theorem. 
The strict inequality > u Q follows from the proof of lcmma l3.3l If = uq > 1, 
then ^ = u + V — is still positive increasing for r > which is impossible. □ 

Appendix A. Transformation to an autonomous system 

Putting u(r) = e~ 2s (/)(s), V(r) - 1 = er 2s W{s) with s = lnr the system flOJ 
transforms into the autonomous system 

4> + (d - 6) <j> - 2(d - 4)0 = W(j> 

W + (d - 6) W - 2(d - 4)V7 = 4> 2 
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subject to the boundary conditions 

(A.2) lim e- 2s 4>(s) = u e R + , lim e~ 2s W/(s) = -1. 

s — > — oo s — > — oo 

For uq G Q we have the asymptotic behavior 

(A.3) lim e- 2s <f){s) = 0, lim e - 2s W(s) = - 1. 

s — >oo s — >oc 

We can associate an energy of the system given by 

(A.4) 2E = 2 - 2(d - 4)0 2 + ^W 2 - (d - 4)W 2 - W(f 

which satisfies 

(A.5) E = -(d-6)(j) 2 + ^W 2 ). 

We should note that though the new system is invariant under translations in s, the 
boundary conditions (|A.2|) break this symmetry and therefore the solutions are not 
translation invariant. However, e~ 2s (j)(s) and e~ 2s W(s) are translation invariant 
which corresponds to the dilation invariance of the original system (|1.7I) . 

If d — 6, then system (| A. 1[) is Hamiltonian and solutions satisfying the boundary 
conditions (|A.2|) have zero energy. We look for a solution such that <fi = — W (i.e. 
uq = 1). Then the zero energy condition reads 

(A.6) <p 2 - Act) 2 + \<t? = 0. 

Since <j> vanishes at ±oo there is so S K such that <j>(so) = 0. Equation (|A.6[) yields 
then 4>(sq) = 6. Integrating (|A.6[) we get 

(A.7) 2(s-s ) = 

J6 



and we obtain the solution 

6 



(A.8) 0(a) = — -j- -. 

cosh (s — so) 

The boundary condition (|A.2[) yields e 2s ° = 24 which after changing variables gives 
the solution ifTTTO]) . 

Appendix B. A Sobolev inequality 

In this part we make the simple observation that the Sobolev-inequality (|1.12[) 
and its optimizers can be obtained by combining two inequalities of [8] . We denote 
by i? 1 (]R !i ) the space of functions / : M. C for which / and V/ are square- 
integrable. 

Proposition B.l. Let d > 4. For all f,g£ H 1 (R d ) 



(B.l) ff \x-y\- A \f{x)\ 2 \g{y)\ 2 dxdy <C d [ \Vf{x)\ 2 dxj \Vg{x)\ 2 dx 
with 

(B.2) C d 



with 

4(d-l) 
d?{d~2) 2 {d-A) 

and optimizers (up to translations and dilations) given by 
(B.3) f{x)=g{x) = {l + \x\ 2 )-^. 
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Proof. We apply inequality (1.1) of [8] with p = t = to get 

// \x- y\- 4 \m\ 2 \9(v)\ 2 dxdy < N pAd \\f\\l p \\g\\% 

and Np^d given by eq. (3.2) of 8 and then the Sobolev inequality (1.2) with 
Kd given in eq. (4.11) of [8]. Both inequalities have the same optimizers which 
concludes the proof. □ 
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